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ABSTRACT 



I present a detailed SUSY QCD calculation of the decay rate of the lightest Higgs 
boson H2 into two gluons, where all quarks and scalar quarks are taken within the 
relevant loop diagrams. I include the mixing of all the scalar partners of the left and 
right handed quarks and show that their contribution is more than several tens of 
per cent compared to the quark contribution in the MSSM for some SUSY parameter 
space. Furthermore the MSSM fermionic contribution is enhanced by several factors 
for large tan (3 and large Higgs masses. As a result, the two gluon decay rate of H2 
is much larger than the two gluon decay rate of an equal mass standard model Higgs 
boson. I further compare the decay mode of H2 —* gg to the similar decay modes of 
H2 cc and H2 —>■ bb including one loop QCD corrections and show that in some 
cases T{H2 gg) is even higher than T{H2 cc), although still much smaller than 
T{Hl-.bb). 
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I. INTRODUCTION 



The Higgs boson is the last particle in the standard model (SM), which yet lacks any 
experimental evidence. Its discovery therefore is of great importance. The instru- 
ments of discovery will be LEP if the Higgs mass is smaller than the Z boson mass 
and LHC for higher masses. While for a Higgs mass smaller than twice of the gauge 
boson mass the most important decay modes for its discovery will he H ^ qq (here 
q — c,b) and ii/" — 77 and to some extent H — > gg, it will be the decay into two W 
or Z bosons for higher masses of the Higgs boson. 

It is well known that the SM is not a sufficiant model when considering unifi- 
cation theories. The favourite model beyond the SM is its minimal supersymmetric 
extension (MSSM) [1] . The content of Higgs particles in the MSSM is quite different 
than the one of the SM: it contains two scalar Higgs bosons Hi,H2, one pseudo- 
scalar ifg and one charged scalar H^. The most important point is that the mass 
of the lightest Higgs particle mjjo has to be smaller than the Z boson mass at tree 
level and is enhanced to a maximum value of around 130 GeV when loop corrections 
are included [2] , thus making a SUSY Higgs boson more experimentally reachable. 

For a Higgs boson far less massive than the Z boson, the H2 — > 77 is the most 
important decay mode and was analyzed in [3] (for the SM) and in [4] (for the 
MSSM). For values of the Higgs mass up to twice the W and Z boson masses the 
decays H2 qq and H2 gg will become more important. The QCD corrections 
to the first decay mode within the SM were considered in [5] (and references therein) 
and within the MSSM in [6]. The second decay mode was considered in [7] and two 
loop QCD corrections within the SM were considered in [8] and found to be rela- 
tively large: about 60%. In this paper I show that the MSSM leads to a much higher 
Higgs into two gluons decay rate than the SM for some supersymmetric parameters, 
making this decay mode more interesting. 

It will be difficult to measure T{H2 — > gg) due to QCD jet background, although 
it might be experimental measurable at future e~^e~- colliders [9]. Therefore it is 
important to consider all kind of models in regards to this decay mode and as I will 
show the scalar quarks contribution can be several tens of per cent compared to the 
quark contribution and also to the H2 — > cc decay rate for some SUSY parameter 
space after summing over all scalar quarks. Furthermore the results presented here 
can also be used for the Higgs production via gluon fusion in pp colliders [7,8]. 

In the next section I present the calculations and discuss the results in the third 
section. In the calculation I include the mixing of all scalar partners of the left and 
right handed quarks, which is expected to be of importance in the top quark sector 
due to the high top quark mass of 174 GcV reported by the CDF group [10]. In the 
bottom quark sector I also include one loop effects. As a surprise I also find, that 
the mixing is not negligible in the charm and strange quark sector independant of 
the value for tan/9, the ratio of the Higgs vacuum expectation values (vev's). 
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II. SUSY QCD CORRECTIONS TO gg 

In the MSSM there are strong relations among the masses and mixing angles of 
the Higgs bosons. Given two values e.g. tan /3 and the light Higgs boson mass m^o 
all the other masses and angles can be obtained by the following equations including 
loop corrections [11]: 
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Here is the soft supersymmetry breaking scalar mass. In eq.(l) I have neglected 
the contribution coming from the bilinear Higgs mass term and A / the parameters 
descriping the strengths of nonsupersymmetric trilinear scalar interactions. Their 
inclusion only changes the results by a few GeV as stated in [11]. 

In the SM the decay mode of the Higgs boson into two gluons occur via one loop 
diagrams with all quarks taken within the loop as shown in Fig.l. After summation 
over both directions of the interior fermions momenta or, which is topological equiv- 
alent, adding the diagram where pi <-> p2 the amplitude is finite and the result is 
given by: 
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with Aq = {'iTiq/mHoY . Before I present the results of the scalar quarks contribution 
to the Higgs decay into two gluons I first want to comment on their mass matrices 
in the MSSM. The mixing term of the scalar partners of the left and right handed 
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quarks is proportional to the quark masses and hence was neglected before the top 

quark was discovered as very heavy. In the calculation here I include the mixing of 
all scalar quarks and present the result in their mass eigenstates, that is instead of 
the current eigenstates qL,R I work with the mass eigenstates 

qi = cosQql + sin Qqji q2 = — sin Qqi, + cos QqR (4) 

In the scalar up quark sector we have the following matrix [12]: 



" V + cot P) ml^ +ml + 0.15Dz J 

Here Dz = cos 2/3, 0.35 — — Bu sin Q\y and 0.15 = Cu sin &w t^ul r soft 
SUSY breaking mass terms, A,^ the parameter from the trilinear scalar interaction 
and 11 the mixing mass term of the Higgs bosons. Here u stands for all three families 
up, charm and top. 

For the mass matrix of the bottom quark, we have to be more careful since it 
is well known that charged Higgsino exchange in left handed scalar down quark self 
mass diagrams leads to a scalar down quark squared mass with a term proportional to 
the up quark mass squared and to flavour changing gluino-scalar down quark-down 
quark couplings [13]. For the mass matrix of the scalar down quark we therefore 
have to take at 1 loop level: 

2 fm^~^+ml-0A2Dz -\c\ml -md{Ad + n tan P) \ 
{ -md{Ad + Attan/3) + - O.OSDz J 

with T3 — Cd sin ©VK = —0.42 and sin Qw — —0.08. Here d stands for all three 
families down, strange and bottom. The value of c is negative and of order 1 (|c] 
increases with the soft SUSY breaking mass term ms and decreases with the top 
quark mass [14]). In the following I take c = —1, although I keep in mind that 
it is more likely smaller. In the calculation it turns out that the mixing of the 
flrst generation is negligible as expected, whereas in the second generation sinQq ~ 
0.1 — 0.5 (the last value only for tan/3 » 1) and therefore not negligible. In the third 
generation sin ~ 1/ \/2 due to the heavy top quark mass. For the scalar bottom 
quark the mixing angle only becomes that big when tan/3 ^ 1. 

In the MSSM we have to add up the two diagrams shown in Fig. 2. After 
summation the amplitude is flnite and the result is given by: 
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with = sin^ Qw and again = m^o on mass shell. Note that the non diagonal 

terms K^2 ^ ^^^Y contribute when the scalar mass eigenstates differ, which 

mainly is the case for the third generation. The amplitudes in eq.(2) and eq.(7) lead 
to the following decay rate: 

- 99) = , , ■ 2r' 2n ^1^^ - 
Stt^ sm"^ Qw cos2 m^o 

Here I take the strong coupling constant Higgs mass dependant as given in eq.(9) of 
Ref. [5]. If Tq is set to eq.(8) reproduce eq.(2.29) given in [15]. Before I discuss 
the results I want to make some comments. The amplitudes given in eq.(2) and 
eq.(7) can also be used when considering Higgs boson production via gluon fusion. 
However there one also has to include the Z boson production via gluon fusion. This 
becomes important when considering heavy lepton production at hadron colliders 
[16,17]. It is also interesting to see that it is possible to produce Z bosons via gluon 
fusion, whereas the Z decay into two gluons is identical to when the Z boson is on 
mass shell = m|, although its amplitude is not (Yang's theorem [18]). 

The scalar quarks do not contribute to the Z boson decay into two gluons (after 
adding the diagram with pi <-> p2 the result is identical to after Feynman integra- 
tion). When quarks are taken within the loop only the terms with 75 survive after 
integration, which also indicates that because of charge conservation no photons are 
produced by gluon fusion (Furry's theorem). The amplitude is then given by: 

'^'^ = cos©^ (ty ^'^^ - + 2A,/Jz6«^,,pfplQ^e-e;fe^ (9) 

with c^^ = Tsq — CqSy^r and c^^ = —^qS^y. Eq.(9) agrees with eq.(l) in [19]. Eq.(2), 
eq.(7) and eq.(9) reproduce the results presented in [17]. When considering the Z 
decay I obtain after squaring and summing over the polarization states of the Z 



4 



boson a term proportional to q^qy^g^'^ — q^q^ /m\) = on mass shell q^ — m^, 
whereas when the Z boson is produced via gluon fusion and then be treated as a 
propagator I have a term proportional to q^{g^'^ — q^ q^ / m\) / {(f' — m^) = —q^/m^ 
and therefore of importance in the gluon fusion process. 

In the next section I will discuss the results of the lightest supersymmetric Higgs 
boson into two gluons decay rate obtained in eq.(8). 

III. DISCUSSIONS 

As input parameters I take mtop = 174 GeV and for the strong coupling con- 
stant the Higgs mass dependant function as = 127r/[(33 — 2n) ln{m'^jjo / -^qcb)] '^^^^ 
n = 5 and Aqcd = 180 MeV. I use ma^ = rriuj^ = m^^ = mj^ = ms = = A^. 
To see how big the contribution of the scalar quarks compared to the quarks is I 
plot in Fig. 3 the ratio T'^+i/T'^ of the decay rate T{H2 — > gg) as function of the soft 
SUSY breaking scalar mass ms for a fixed value of fi = 250 GeV and two different 
values of the Higgs mass 711^0 — 60 GeV and 120 GeV and 3 different values of 
tan/5 = 3 (solid line), 10 (dashed line) and 60 (dotted line). Higher values of tan/? 
are preferred in superstring inspired Eq and 50(10) models. For F^ and F^+^ I have 
taken Tq as given in eq.(2) with the couplings K'^^^ , that is including the large en- 
hancements (relative to the SM) due to large tan/3. As a result I have that for small 
values of ms the scalar quarks contribute even more than the quarks, although their 
contribution decrease rapidly and remains only a few per cent for ms > 600 GeV. 
For tan/3 = 60 the scalar quarks contribution diminishes the ratio for ms < 350 
GeV and enhances it for higher values. The influence of is very small for small 
tan/3 and becomes more important for very high tan/3 values. For small values of 
tan/3 higher values of /x enhance the decay rate a little bit. For high values of tan/3 
it is the other way around and the differences are larger. A negative value for fj, leads 
to a bit smaller values of the decay rate. 

In Fig.4 I have plotted the ratio of the Higgs into two gluons decay rate of the 
MSSM compared to the SM, that is F^+'^/F'^, where I have taken F"^ as it is in the 
SM that is without the couplings K'^^^ ^ whereas I included them in F*^"*"^. For the 
Higgs mass I have taken 60 GeV. As a result I have that for scalar masses smaller 
than 500 GeV the Higgs into two gluons decay rate is enhanced by several tens of 
per cents in the MSSM and gives the same result than the SM for higher values 
of the scalar masses. As in Fig. 3 for tan/3 = 60 and ms < 350 GeV the decay is 
diminished. 

In Fig. 5 I have done the same as in Fig.4 but for a Higgs mass of m^o = 120 
GeV. Here the results are quite different than in Fig.4. As a result I have that in 
the MSSM the Higgs into two gluon decay rate is enhanced by several tens of per 
cent for tan/3 = 3, by a factor of 2 — 3 for tan/3 = 10 and ms < 300 GeV and even 
by an order of magnitude for tan/3 = 60 with the highest contribution for a scalar 
mass around 550 GeV. 

As I have shown in Fig. 3 the scalar quarks decouples for ms > 600 GeV. The 
reason why the branching ratio as shown in Fig. 5 is still larger than 1 even for higher 
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values of the scalar mass is that is quite different in the MSSM with the couplings 
K'i^^ than it is in the SM without these couplings. In the SM the main contribution 
is from the heavy top quark and a few per cent from the bottom quark. The contri- 
bution of the other quarks are negligible due to their small masses. In the MSSM the 
bottom quark contribution becomes as important as the top quark contribution for 
large tan/3 values eg. the ratio Tlj^/F^ggj^ becomes very small depending on the 
size and sign of sin a, which becomes relative large around ±0.5 and thus leading to 
large values of r'^+'^/F'^ as seen in Fig. 5. For a small Higgs mass of 60 GeV as I have 
taken in Fig. 4 sin a remains always smaller than around —2 x 10~^ and therefore 
keeps the bottom quark mass contribution as small as in the SM. 

Some curves in Fig. 4-5 start at different values of ms because, for values of ms 
higher than 600 GeV I obtain an unphysical negative mass squared for the pseudo 
particle if tan/3 = 3; for tan/3 = 10 the unphysical region is when ms — 650 
GeV; whereas for tan/3 = 60 mjjo is physical for all ms- 

A negative eigenvalue of the scalar bottom quark mass also occurs if ms < 200 
GeV for tan/3 = 3 and 10 or ms < 300 GeV for tan/3 = 60. Here the parameter c in 
the scalar bottom mass squared matrix eq.(6) is of importance, neglecting it would 
allow us to use ms as small as 100 GeV (for tan/3 = 3 ) without running into one 
negative mass eigenvalue of the scalar bottom quark mass, with the result that F^"*"^ 
can become much larger than F^. Unfortunately even for smaller values for c ~ —0.5 
I obtain negative values with such a small scalar mass. Since c cannot be neglected 
when including loop corrections I excluded those regions in the figures. 

In Fig. 4 and Fig. 5 I also have plotted the ratio of the decay rates F(i/2 — 
gg)/T{H2 cc). For the decay rate T{H2 hh,cc) I used the tree result includ- 
ing the SM QCD corrections as given in eq.(8) of Ref. [5] with the changes of the 
tree level couplings within the MSSM. I did not include the SUSY QCD correction, 
because they are far smaller than the SM QCD correction as I have shown in [6]. 
There I showed that for tan (3 — 1 SUSY QCD corrections do not contribute at all to 
this decay mode (sin(Q; -\- j3) — and presented the results in the limit of tan /3 3> 1. 
There I did not include the mixing of the scalar charm and bottom quarks, but since 
mg^ ~ m^.^ even a large mixing angle will not change the results for the decay rate 
r{H2 — > cc) presented there. This might not be true for T(H2 bb) which I did 
not consider here since it is much higher than T{H2 gg), by a factor of at least 
50. Therefore in Fig.4 and Fig. 5 I only compared the Higgs into two gluons decay 
rate with the one to the charm- anti-charm quarks. In Fig.4, the dependence of the 
ratio to the scalar mass ms is very small, since F*^"*"*^ becomes very small and the 
dependence of F^'^ to ms is only via the angles cos^ a/ sin^ /3, which is compensated 
by the K^^^ coupling in the Tg term. For a very large scalar mass the ratio remains 
constant with a value of around 0.31 independant of tan /3. A quite different result 
I obtain in Fig.5, especially again for tan/3 = 60, for the same reason as explained 
above. The shape of the figures is quite similar compared to the ratio F'^+'^/F'^. For 
scalar masses much higher than 1 TeV the ratio remains constant with a value of 
around 1 independant of tan /3. 

Finally in Fig. 6 I show the influence of the Higgs mass to the decay rate 



6 



T{H2 gg) for a fixed value of = 250 GeV and ms — 300 GeV and three 

different values of tan/3 = 3,10 and 60. In the case tan/3 = 31 obtain negative 
values for the mass squared of the pseudo Higgs in the range of 95 < m^o < 105 

GeV which therefore has to be excluded. As a result I have that r^+^/F'^' is weakly 
dependant of the Higgs mass, not so r^+'^/F'^'^, which shape is basically dominated 
by cos a and sin (3. 

IV. CONCLUSION 

In this paper I presented the corrections to the lightest MSSM Higgs boson de- 
cay into two gluons when scalar quarks are taken within the loop. I included in my 
calculation the mixing of all scalar quarks although it only becomes important for 
the second and third generation. I have shown that scalar quarks lead to a decay 
rate of the same order as the quarks in the SM for values of ms smaller than 600 
GeV. In the SM the largest contribution comes from the top quark due to the 
in Tq. In the MSSM the Tq are of the same order for all scalar quarks and therefore 
contribute many more terms to T{H2 gg) than the SM alone. Furthermore in 
the MSSM the Tq can become much larger than in the SM for tan /3 » 1 and large 
negative or positive sin a. I also have shown that the Higgs into two gluon decay rate 
can become even larger than the decay into charm- anti-charm quarks for tan (3 — 3 
and the Higgs mass larger than around 80 GeV and for tan/3 = 10 and 60 and the 
Higgs mass larger than the Z boson mass, but still remains more than a factor of 50 
smaller than its decay into bottom- anti-bottom quarks. 

Although the decay of the Higgs into two gluons will be difficult to measure 
it is of importance to know how big the influence of models beyond the SM might 
be. Furthermore it might be measurable at future e~^e~- colliders [9]. The ampli- 
tudes given here can also be used when considering Higgs and Z boson production 
in hadron colliders via gluon fusion. 
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FIGURE CAPTIONS 

Fig.l The penguin diagram with up and down quarks within the loop leading to the 
H2 99 decay. 

Fig. 2 The penguin diagrams with scalar up and scalar down quarks within the loop 
leading to the H2 — * 99 decay. 

Fig. 3 The ratio r^+i /T'^ of the decay H2 gg a,s a function of the soft SUSY breaking 
scalar mass term ms for a fixed value of fi = 250 GeV and rrijjo = 60 GeV and 
120 GeV and three different values of tan(3 = 3 (solid line), 10 (dashed line) and 
60 (dotted line). For tan/? = 3 the upper curve is for 171^0 = 120 GeV and the 
lower curve for 60 GeV. For tan (3 = 10 and 60 it is the other way around. Here 
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I have taken for the equation as it is given in eq.(2) including the couphngs 

Fig.4 The ratio r^+'?/r'? and T^+^/F'^^ of the decay ^ gg as a function of the 
soft SUSY breaking scalar mass term ms for a fixed value of n — 250 GeV and 
uifjo = 60 GeV. tan/? as in Fig. 3. Here I have taken for F^ the equation as it is 

given in the SM, that is without the couplings K'^^° . 
Fig.5 The same as Fig.4 with mfjo = 120 GeV. The upper curves are for T^+i/T'^'^. 
For tan/3 = 3 and nis > 600 GeV the pseudo Higgs becomes a negative mass 
squared. For tan/3 = 10 the same happens for a small region when ms ~ 650 
GeV. 

Fig. 6 The same as Fig.4 but as a function of of m^o for fixed values ms — 300 GeV 
and iJ, = 250 GeV and tan/3 as in Fig.3. Here the upper curves at the higher 
Higgs masses are for T^+i/T'^^. In the case tan/3 — 3 the pseudo scalar mass 
squared m'Lo becomes negative in the range 95 < mjjo < 105 GeV. 
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